Resistivity due to low-symmetrical defects in metals 
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! Abstract 

CN ' The impurity resistivity, also known as the residual resistivity, is calculated 

ab initio using multiple-scattering theory. The mean-free path is calculated 

^ ■ 

in 

qq ' by solving the Boltzmann equation iteratively. The resistivity due to low- 

^1 ' symmetrical defects is calculated for the FCC host metals Al and Ag and the 

00 

o> 



BCC transition metal V. Commonly, j noise is attributed to the motion of 



& ' such defects in a diffusion process. The results for single impurities compare 

a" 

well to calculations by other authors and to experimental values. 



I. INTRODUCTION 

The theoretical explanation for the electrical resistivity is well-known. Electrons move 
through a regular lattice of metal atoms without any resistance. As soon as irregularities are 
introduced in this metal electrons are scattered, which gives rise to a finite resistivity. The 
temperature dependence of this quantity is mainly due to scattering of electrons by phonons. 
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At zero temperature, when no phonons are present, the resistivity is determined by defects 
only, such as impurity atoms. Then it is the only remaining contribution and therefore it 
is often called the residual resistivity. In this paper the resistivity due to impurity atoms 
embedded in the metal lattice is considered, the impurity resistivity, which is extensively 
studied experiment allyS. 

An interesting problem is the problem of resistance noisei. Over a large range of frequen- 
cies the spectral density varies as 1/ /. This can be explained, if these resistance fluctuations 
arise from a kind of diffusion process. In most cases the frequencies range from 1 to 1000 
Hz, which correspond to typical times between jumps. The noise is attributed to a defect, 
which can be of any kind, jumping back and forth. A simple example of such a defect 
is an impurity-vacancy pair, of which we are able to calculate the resistivity for different 
orientations. 

A lot of attempts have been made to calculate the impurity resistivity. The simplest 
methods consider an atom or a cluster of atoms embedded in free spaceiS. More sophisti- 
cated approaches use ab initio methods like the Korringa-Kohn-Rostoker theory (KKR)H 
to describe an impurity embedded in a metal lattice. If this formalism is applied for two 
spin directions magnetic impurities and materials can also be treatedi. In most 
substitutional or interstitial&i impurity atom is considered. In this work we mainly concen- 
trate on the resistivity due to defects, playing a role in substitutional electromigration, such 
as a vacancy, an impurity-vacancy pair and an atom on its way to a neighbouring vacant 
lattice site. The symmetry of most of the considered defects is reduced compared to a single 
impurity atom, which magnifies the required computational effort. 

The theory, which is used to calculate the impurity resistivity is described in Sec. ||. 
The theory makes use of the calculation of the electron wave function described by Dekker 
et aiB, which already requires a heavy computation of a Green's function matrix. In 
Sec. |l Jl| results are shown for the host metal Al. The calculations for single 3d and Asp 
impurities in Al are compared with experimental and other theoretical values in Sec. A| . 
In Sec. 1 1 [B] various calculations are reported, which are interesting in view of the reliability 



measurements mentioned above. Vacancies and moving host atoms in Al are considered. 
Resistivity calculations for impurities, a vacancy, several impurity- vacancy pairs and and an 
impurity at the saddle point in the FCC metal Ag are done in Sec. [TV]. Results from similar 
calculations for the BCC transition metal V are reported in Sec. |V|. A summary is given in 
Sec. jV]. 



II. THEORY 

First the general theory will be presented. After that some equations are given for the 
resistivity due to low-symmetrical defects. Finally a new expression for the generalized 
Friedel sum, used in the present paper, will be given. The conductivity of a sample can be 
calculated performing an integration over the Fermi surfaced 

in which the velocity of an electron with quantum numbers k = (nk) is extracted from 
the host electronic structure. A finite electron mean free path is due to the presence of 
defects or phonons and can be calculated by solving the equation 



A, = rl 



v fc + J! Afc'-Pfc'fe 
k' 



(2) 



This equation follows easily from the linearized Boltzmann equation. In this paper scattering 
by a static defect is considered. The defect can consist of a number of perturbed host atoms, 
an impurity and one or two vacancies. The probability rate P^k for the transition through 
scattering from state k to k' determines the electron lifetime r° 
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k 
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For a low concentration c of a certain kind of defect, the transition probability P^k for elastic 
scattering is given by 



P fc / fc = 2>KcN\Ty k \ 2 5(e k - e k >). (4) 



The calculation of the transition matrix T^k requires knowledge of the electronic wave 
function of the alloy. This wave function can be calculated using multiple-scattering theory. 
The formulation of this theory is given by Dekker et a/.0. For the sake of clarity, some 
quantities appearing in the theory, which are necessary in the evaluation of the impurity 
resistivity, will be given here too. 

The alloy wave function coefficients c^l and host wave function coefficients c\ nL are 
related by a matrix equation, 

CknL = E ^LL' C kn'L'- (5) 
n'L' 

The matrix label n refers to an atomic site, either at a host position Rj or at an alloy 
position Rp, and L = (l,m) summarizes the angular momentum labels. The matrix A 1 ^, 
will be defined below. The host wave function coefficients are evaluated at the Fermi energy 
E-p = k, 2 and can be written as 

^° L (k)e» kR " 

knL ^(-(d\ /de k ))^ [ } 

The vector W^ L (k) is defined by 

JW^M = £ hu(K Rn)/V^(k), (7) 

L' 

where 6(k, R n ) is a lattice sum 

b(k,R n ) = J2B(R nf )e-^> (8) 

f 

and i l Vl and X q are an eigenvector and the corresponding eigenvalue of the KKR matrix 
M (k) = t h 1 — 6(k, 0). The matrix B is defined with Gaunt coefficients Cll'l" and spherical 
Hankel functions h~^(r) = hf{Kr)Yi{r) as 

B LV {R) = i l - l '- l KY,C L L'L"fh + L {R). (9) 

L" 

It has to be stressed that the lattice sum in Eq. (|5p extends over all host positions when 
R n is not a host lattice position. When it is a host lattice position Rj, the corresponding 



term is excluded. The q = label in Eq. fl6|) refers to the eigenvalue, which corresponds to 
a zero KKR matrix and thereby determines the electronic structure of the metal. 
The matrix Aj^i in Eq. ([5]) is defined a 

Alt = £ - G v ° idt )~ in Z i 1 ~ g ™ dth )lZ> (10) 

niL\ 

where the scattering matrices of the atomic host potentials t 1 ^ and the ones of the atomic 
alloy potentials t n are calculated from their phase shifts 

tl = -sin(r] nl y^. (11) 

The host phase shifts for an alloy position p, rj^, are defined to be zero, if the position p 
does not coincide with a host position. The alloy phase shifts for the host position j are 
defined to be zero if the position does not coincide with an alloy position. 

The formalism is made suitable to handle more general defects by making use of a void 
system as a reference system instead of the unperturbed host. The impurities and perturbed 
host atoms are replaced by free space in this reference system. The Green's function matrix 
of this reference system is calculated from the host Green's function matrix 

£void,n„' = gnu' _ £ gnh^ + Q)^Q^' . (12) 

The host Green's function matrix is calculated by an integration over the Brillouin zone 

G nn> = ^— I rf 3 fc[6(k,R„ n 0+Kk,Rn)M- 1 (k)6 T (-k,R n ,)]e 4kR ""'. (13) 

As derived by Van Ek et a/.i and Lodder et aZ.0 T^k can be written within multiple- 
scattering theory as 

T k 'k = E C fc'nL T L C knL, (14) 
nL 

where T£ is defined as 

T n L = --sin( Vnl - Vni)^ 1 '^- (15) 

K 
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We define a new quantity QknL as 

QknL = —F=-T1CknL- (16) 



I- 1 



K 

Now the sum over k' in Eq. @ can be rewritten as a Fermi surface integral and a set of 
equations in terms of W® L (k) and QknL can be derived straightforwardly. The equation for 
Afc becomes 



A fc = r° k 



Vfc + C ^2 QknLQkn'L'^LL' 



(17) 



nn'LL' 

where I is a Fermi surface integral with as a factor in the integrand 



w 2tt /- ^ z (k)*A fc < a ,(k) & . R ^ ; 
^ B z JFS |VfcA u (k)| 

Eq. (|17D can be solved iteratively. In the calculation of r° we can make use of the optical 
theorem, which states that the sum over k' in Eq. ([3D can be connected to the diagonal 
element of the transition matrix 

7-r 1 = -2clmT kk . (19) 

The comparison of the two expressions for r°, (§) and (p^9[ ), can serve as a test for the 
accuracy of the Fermi surface integrals. For a more complete description of the theory for 
host and alloy wavefunctions, the reader is referred to Dekker et a/.0. Here we just add, that 
an initial A has to be inserted in Eq. (|18D , e.g. Ak = t^v^ or the Ziman approximationS. 
This leads to a new set of A^ according to Eq. ([17]) . With this new set the integrals in Eq. 
(|18|) can be recalculated. This procedure is repeated until the new set equals the inserted 
set. 

Now we give the current density-field relation for a metal containing low-symmetrical 
defects. In such a metal the resistivity is anisotropic, i.e. it depends on the direction of 
the current. So, the relation between the electric field and the current density for e.g. an 
impurity-vacancy pair in the FCC structure is given by 

j = ^E,| + + — E z , (20) 

p\\ p± pz 



where En lies along the jump direction of the migrating atom and both Ej_ and E z are in 
a perpendicular direction. The different directions are shown in Fig. For an impurity- 
vacancy pair in the BCC structure there are two inequivalent directions, which are displayed 
in Fig. ^| and therefore the current density can be written as 

j = -E,| + — E ± . (21) 

Eqs. (|20|) and ([H]) describe the current density in a sample, containing only one kind of 



defect, with one particular orientation. In a real sample the orientations of a defect are 
distributed randomly. Such a distribution results in a scalar resistivity, which is given by 

PFcc = U- + — + -) (22) 



3 \p\\ P± Pz 



for an FCC metal and by 



Pvhc = \(± + f ) (23) 
3 \P\\ P±J 

for a BCC metal. 

Finally, in order to check the requirement of charge neutrality for the potentials to be 
used, we need an expression for the generalized Friedel sum. We will show that it is possible 
to derive such an expression, using the formalism presented above. According to Lodder and 
Braspenning@ the electron density of states of a system n(E) can be written with respect 
to an arbitrary reference system as 

n(E) = n ief (E) + -Im— Tr In T(E), (24) 
7T dE 

where T(E) is the t matrix of the system, with respect to the reference system. Convention- 
ally the unperturbed host has served as a reference system for a dilute alloy. For a general 
defect the void system serves as the natural reference system. In that case the t matrix of 
the system can be written as 

T(E) =t{\ -g void t)-\ (25) 
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The integrated density of states N(Ep) = f F n(E)dE up to the Fermi energy counts the 
total number of electrons accomodated in the system. The difference in the number of 
electrons between the alloy and the host, Zp, is found by subtracting N host (Ep) 

Z F = N{E F ) - N host {E F ) = 
2 2 

-arg det t - -arg det(l - G vold t) + 

7T 7T 

--arg det t h + -arg det (1 - £ void t h ), (26) 

7T 7T 

which is the generalized Friedel sum. In the case of spherically symmetric scatterers this 
general expression simplifies to 

z F = - 52(21 + i)n? - -iv c h lustcr 52(21 + 
n pi n i 

--argdet(l - £ void i) + -arg det (1 - £ void t h ), (27) 

TT 7T 

in which A' r ( ! 1 luster is the number of host atoms in the void region. 

This expression is more general than expressions used in the past,0 which only apply to 
simple substitutional and interstitial alloys for which no intermediate void reference system 
was needed. We will show that Eq. (p^) reduces to well-established expressions applicable 
to those simple systems. In order to do this, it is useful to extend the sum in Eq. (|T2D to 
interstitial sites. This can be done by defining host scattering matrices for those positions 
as = 0. By that the elements of the matrix (t h 1 + Q) = t h (l + £t h ) -1 are equal to 
zero, when one of the two or both indices refer to an interstitial site. The resulting matrix 
equation g vold = Q — g(t h 1 + Q) Q contains only matrices of the same dimension, and Eq. 
(O) can be rewritten as 

g = (i - ^ void t h ) _1 ^ void . (28) 

Note that this equation can be derived directly from Eq. (|T2| ) in the case of a substitutional 
alloy, where only lattice sites are occupied. In the case of an interstitial impurity the matrices 
are enlarged due to the presence of the interstitial atom. 
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The addition of a non-scattering atom does not affect the host charge. This can be seen 
from Eq. (|2T|), and is trivial from a physical point of view. The matrices of the third and 
fourth term can be multiplied, leading to 

(1 - g void t h y\l - g void t) = 1 - (1 - g™ id t h )~ l g void (t - t h ) = 1 - Qit - t h ). (29) 

Hence, the Friedel sum is given by 

Z F = - "£(21 + l)( V f - rjt*) - -argdet(l - Q{t - t h )), (30) 

pi 

which has been applied in the past to substitutional^ and interstitial^ alloys. 



III. IMPURITY RESISTIVITIES IN AL 

A. 3d and 4sp impurities in Al 

In this section a single 3d or Asp impurity is considered embedded in unperturbed Al 
host. This means that the charge transfer to the surrounding host atoms as well as lat- 
tice distortion are neglected. Furthermore, an impurity atom has an assumed electronic 
configuration, which in reality may depend on its metallic environment. From Fig. |3|, in 
which the calculated impurity resistivities are shown, it is clear that this configuration is 
very important. The filled circles refer to calculations in which the impurity atom has one 
4s electron. The values indicated by filled squares are obtained for impurity atoms with 
two 4s electrons. The impurity resisitivity of atoms having two 4s electrons decreases with 
increasing atomic number, while it shows a maximum for Mn, when only one 4s electron 
is present. The experimental values^, indicated in the figure by asterisks, also show such a 
maximum, but the values are underestimated by the calculations. 

The potentials used in the calculations just described do not lead to a charge neutral 
system, which is unphysical. The neutrality can be restored by adding a surface charge to 
the atomic spheres.0 This procedure is called the shifting procedure, because it corresponds 
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to a shift of the atomic potential by a constant energy. The charge of the system is cal- 
culated using the generalized Friedel sum expression given in Sec. [TI|. This procedure has 
been applied to the transition metal impurities with the 3d n 4s 1 electronic configuration, the 
Ca(4s 2 ) and the 4sp impurities. The impurity resistivities, obtained with these potentials, 
are given by open circles in Fig. |3|. The addition of charge leads to an increase of the resis- 
tivity in all cases, except for Sc, Ge and As. The agreement with the experimental values 
becomes much better. For all 4sp impurities and for the transition metal impurities with 
more than six 3d electrons the agreement is very good. 

The addition of surface charge is a crude attempt to simulate the effect of charge re- 
laxation in the alloy. Still, in the case of the 3d impurities Fe, Co and Ni it enhances the 
accuracy of the resistivity significantly. Unfortunately, this is not the case for the other 3d 
impurities. Apparently, the surface charge does not simulate all effects of charge relaxation 
in the right way. Therefore it would be very interesting to repeat the calculations for Sc, 
Ti, V, Cr and Mn with self-consistently calculated potentials. The method of calculation of 
the resistivity is not affected by the use of such potentials. 

The resistivities of these impurities in Al were already calculated by Boerrigter et a/.i, 
Schopke and MrosanB and recently by Papanikolaou et al%. Schopke and MrosanEl used 
the spherical band approximation, which means that the Fermi surface is approximated by 
a sphere. They found resistivities, which were approximately equal to the ones following 
from the well-known free-electron formula of FriedelEl, which only contains the scattering 
phase shifts. Just as the other authors mentioned they found an underestimation of the 
resistivities, which was attributed to the anisotropy of the Fermi surface. Papanikolaou et 
al% tried to incorporate these anisotropy effects in a tricky way and found values for the 3c? 
impurities, which were too large. In our calculation this anisotropy is fully and consistently 
taken into account, but still the impurity resistivities are underestimated. 
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B. A migrating Al atom 



According to our calculation the resistivity of a vacancy in Al is 0.57 /if2cm/at%. We 
used host phase shifts for all surrounding Al atoms. In first order the resistivity is the sum of 
the resistivities of the separate scatterers. Therefore it is likely that the vacancy resistivity is 
underestimated. In the present case account of the scattering by the first shell enlarges the 
resistivity only slightly, to 0.60 /xf2cm/at%. Our value contradicts with earlier calculations 
of Van Ek et who found 0.93 /if2cm/at%. 

The vacancy resistivity is also extracted from simultaneous measurements of the resistiv- 
ity and the expansion of both the total volume and the lattice constant in an Al sample^. 
In this way a value of 3.0 //f2cm/at% is found, which is much larger than the value we found. 
This could have several reasons. One of the reasons can be that the electronic structure of 
the vacancy defect is not calculated self-consistently. From the previous subsection indeed a 
strong dependence on the electronic structure was observed. Another reason may be that the 
volume expansion is not entirely due to the absorption of vacancies or that the enlargement 
of the resistivity is not merely due to the presence of vacancies. 

During a jump the resistivity changes from the initial value, via the value at the saddle 
point, back to the initial value. The saddle point value depends also on the direction of the 
jump with respect to the direction of the current. In the calculation a single saddle point 
atom is taken into account, so scattering by the two small moon-shaped vacancies next to 
the atom is neglected. This procedure leads to a resistivity which is smaller than the one 
of the vacancy for all directions of the current, namely pn = 0.55/ificm/at% and p± and 
p z both have the value of 0.36/if2cm/at%. The resistivities for the different directions are 
defined by Eq. It is expected that the small vacancies contribute considerably to the 
resistivity, leading to a value, which is larger than the vacancy resistivity. 

Calculations for a pair of vacancies show that the resistivity, averaged over all current 
directions, is equal to the resistivity of two single vacancies. Perhaps a larger cluster of 
perturbed host atoms or self-consistently calculated phase shifts could alter this conclusion. 
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The symmetry of a pair is the same as the symmetry of an atom at the saddle point. 
Therefore Eq. ( pop holds. The parallel resistivity p» turns out to be 0.94 /xf2cm/at%, which 
is considerably smaller than the resistivity in the other two directions (p± = 1.24/if2cm/at% 
and p z = 1.31yuf2cm/at%). The much smaller resistivity of a pair of vacancies aligned along 
the current is easily explained intuitively with the help of Fig. [|. Assuming a monotonic 
relation between the geometrical and scattering cross-sections, the scattering cross-section is 
obviously larger when the pair of vacancies is aligned perpendicular to the current. However, 
from the results for impurity-vacancy pairs, to be presented below, it follows that this 
intuitive, classical explanation does no justice to the quant ummechanical character of the 
scattering process. Microscopically, one has to consider the scattering probability due to 
a pair of potentials v and w, lying at a distance R, which, of course, is not simply equal 
to the sum of the individual probabilities too. Even in lowest order in the potential, this 
probability Pk>k, calculated in the free electron model, so using plane waves, is proportional 
to 

Pwk ~ v\, k + w\, k + 2v Wk w Wk cos((fc' - fc).R), (31) 

in which v k >k — 4vr / r 2 drj Q (\k' — k\r)v(r) is a real quantity for a spherical potential in free 
space. For a pair of vacancies v = w. It is clear that the cosine term does not have a definite 
sign, and that the contribution will be different for different alignments of R. Our results 
for the pair of vacancies imply, that the average contribution of this term is positive for R 
perpendicular to the current, and negative for alignment along the current. For large values 
of R this term will average out, and the individual probabilities just add. 

IV. 5SP IMPURITIES IN AG 

The experimentally obtained resistivities of the 5sp impurities in AgH have already been 
used in the analysis of their wind valence by Dekker et aZ.0. In this section the impurity 
resistivities will be calculated for a single impurity, one next to a vacancy and one at the 
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saddle point during a diffusion jump. In most of the calculations the perturbation of the 
surrounding host atoms is not taken into account. In Fig. [5] it is seen that the calculations, 
indicated by filled circles, and the measurements, indicated by asterisks, show the same 
trend. However the measured values are larger. Only the value of 1.18 pQcm/ai% for the 
4d 10 impurity Pd is an overestimation. A much lower value of 0.02 pQcm/at% is found, 
when a 4c? 9 5s 1 electronic configuration is used for the Pd atom. The experimental value of 
0.44 /zffcm/at% lies between the two theoretical values, which suggests that the electronic 
configuration is a mixture of both. The calculated resistivities are only slightly affected by 
taking into account a shell of perturbed host atoms. A maximum increase of 0.04 /if2cm/at% 
is found for In. 

The shifting procedure to achieve charge neutrality is also applied in this case. Missing 
charge had to be added to the impurity. The resulting values are indicated by open circles 
in Fig. [5[ Just like in the case of impurities in Al the resistivities are enlarged. However, 
the agreement with experiment does not improve in this case, because the enlargement is 
too strong. 

Similar calculations have been done by Vojta et a/.i using self-consistent single-site po- 
tentials. Their results are comparable to ours, but they agree somewhat better with the 
experimental values. This could be the result of a larger muffin-tin radius, they used. Our 
muffin-tin radius is bounded, because of the decreased space at the saddle point. Neverthe- 
less our values are reasonable. 

The resistivities for 5sp impurity-vacancy pairs are given in Fig. || The resistivity for a 
single vacancy is 0.82/if2cm/at%, which is the value for Ag in the figure. The resistivity of 
an impurity- vacancy pair, being aligned with the current, p», is larger than the resistivity, 
when they are aligned perpendicular to the current, p± and p z . This is in contradiction 
with the intuitive explanation for the resistivity of a vacancy pair in Al in the different 
directions in terms of a geometrical cross-section, which is given in Sec. [Ill B| and illustrated 



in Fig. ^. However, this behavior can be understood from the simple expression fl3"Tp. The 
impurity potential w is certainly attractive, which corresponds to an overall negative sign, 
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and a vacancy potential v is repulsive. So, on the average, the cosine term in Eq. (|3~1~D has 
the opposite sign compared with the scattering by two vacancies. This implies a conversion 
of the behavior, in agreement with or finding for the impurity-vacancy pair. Notice also, 
that the resistivity of an impurity- vacancy pair, averaged over all current directions, ^average, 
does not equal the sum of the separate resistivities of vacancy and impurity. The latter sum 
rather equals p\\. 

In Fig. [7| the impurity resistivities at the saddle point for the different current directions 
are compared with the corresponding resistivities for the impurity-vacancy pair. The saddle 
point resistivity follows roughly the one at the initial position. Again p\\ is the largest, but 
for an atom at the saddle point the cross-section is not expected to depend strongly on the 
direction, because the current sees one scattering atom from all directions. Just like in the 
case of Al, the two small moon-shaped vacancies around the saddle point atom are not taken 
into account, which is expected to lead to an underestimation of the resistivity. 

V. TRANSITION METAL IMPURITIES IN V 

The measured resistivities of the 3c? impurities Ti and and the calculated ones of 
Sc, Ti, Cr and Mn in V are given in Fig. |§. The calculated values are lower than the 
experimental values, although the value for Cr lies fairly close. The Mn resistivity is much 
larger than the other ones. The value measured for the hd impurity Ta of 1.5/xficm/at% is 
very close to the calculated value of 1.3/if2cm/at%. 

The calculated resistivity of a vacancy in V is larger than of any of the 3d impurities, 
namely 4.94/if2cm/at%. This results in resistivities of impurity- vacancy pairs, varying from 
5 to 9 pQcm/at%, as can be seen from Fig. |9|. The large value for the Mn impurity is also 
seen in the 3d series in the left panel of the figure, but the effect is not as pronounced as in 
the case of a single impurity. The resistivity turns out to be fairly isotropic, i.e. p\\ ~ p± in 
Eq. 0. 

It is seen that the resistivity of a 4d impurity next to a vacancy tends to be larger than 
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that of a 3c? impurity and smaller than that of a 5c? impurity. The resistivity for the 3c? 
impurities is the lowest for V, while for the 4c? impurities it is lowest for Mo, which has an 
additional valence electron. For the 5c? impurities the resistivity of the impurity-vacancy 
pair decreases monotonically with the atomic number. 

The resistivities for impurities at the saddle point are depicted in Fig. |TU[ They show 
a larger anisotropy. Exceptions are Cr, Mo and W. Apart from the high value of Cr, the 
resistivity seems to decrease monotonically in all three series. The low value for Mn is 
striking in view of the high values for the single impurity and the impurity- vacancy pair. 
The saddle point resistivities are larger than the initial point values. The small vacancies 
on either side of the atom could even enhance this effect. 



VI. SUMMARY 

In this paper a multiple-scattering method has been described for the calculation of the 
impurity resistivity. It makes use of the calculated wave function coefficients, introduced by 
Dekker et The linearized Boltzmann equation can be solved iteratively. One iteration 
step involves the calculation of a Fermi surface integral. The integrand is the product of the 
vector mean free path, which depends on the crystal momentum, and two host wave function 
coefficients. In its present formulation, the method is suitable to handle complicated defects 
such as an atom during a diffusion jump. It has been used to calculate the resistivity due 
to impurities, vacancies and pair defects in Al, Ag and V. 

The resistivities of 3c? and 4sp impurities in Al have been calculated, basically in order to 
see if the calculations make sense. This series of impurities was investigated before by several 
authorHHB and experimental values are available^. Their calculated resistivities turn out 
to depend strongly on the atomic electronic configuration, which is used to construct the 
crystal potential of the alloy. This is especially important for transition metal impurities, 
where e.g. the energies of 3c? and 4s levels are almost equal. In this series it is seen that 
the resistivity decreases with atomic number, when the impurity has two 4s electrons. The 
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shape of the experimentally observed peak is reproduced, when the impurity carries one 4s 
electron. 

Another consequence of the construction of the potentials, the lack of charge neutral- 
ity, can be repaired by adding surface charge to the atomic sphere of the impurity. This 
procedure enlarges most calculated values and improves the agreement with experiments. 
Especially for transition metal atoms with many d electrons, and for the 4sp impurities the 
agreement becomes very good. Apparently that the calculation takes the essential features 
of the scattering process into account. The strong dependence on the electronic configura- 
tion as well as on the addition of surface charge make it interesting to use self-consistent 
potentials in our calculation. 

A vacancy plays an important role in the diffusion process. Its calculated resistivity in Al 
of 0.6 /iOcm/at% is much smaller than the experimentally obtained value of 3 /zf2cm/at%. 
The resistivity of a host Al atom, halfway along its jump path to a neighbouring vacant 
site, depends on the direction of the electrical current and it is different from its value for 
the atom at its initial position. Both the direction and position dependence give rise to 
fluctuations in the resistivity on a timescale of 10 -13 s. The value of 0.41 /ificm/at%, which 
is the average over all current directions, is smaller than the value at the initial position, the 
latter being equal to the resistivity of a vacancy. In this calculation the two small moon- 
shaped vacancies next to the jumping atom are not taken into account and it is expectable 
that they will enlarge the resistivity. The resistivity of a pair of vacancies depends on the 
direction of the current. If the pair is aligned with the current, the resistivity is smallest. 
This can be attributed to a smaller cross-section for such a configuration. If the resistivity 
is averaged over all current directions it equals the resistivity of two single ones. 

The calculated resistivities due to the 5sp impurities in Ag show a similar dependence 
on atomic number as the experimental values.0 Just as for impurities in Al the resitivities 
are underestimated. However, after achieving charge neutrality by adding a surface charge 
to the impurity, they become too large. The resistivity due to an impurity-vacancy pair is 
smaller than the sum of the impurity and vacancy resistivities. When the pair is aligned 
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with the current, the resistivity is largest and approximately equals that sum. The fact that 
the resistivity is largest in that direction is in contradiction with the smaller geometrical 
cross-section. An impurity halfway its jump path has a larger resistivity than the impurity- 
vacancy pair in spite of the neglected small vacancies. 

The calculated resistivities of the impurities Cr and Ta in the BCC transition metal V 
agree fairly well with experiment, while the one of Ti is underestimated. The values for a d 
impurity- vacancy pair and an impurity halfway its jump path are larger than the ones for a 
single impurity. 

In conclusion, it has been shown that the resistivity due to low-symmetrical defects can 
be calculated accurately. The calculated impurity resistivities compare reasonably well with 
the available experimental material. They may even improve when self-consistent potentials 
for the alloy are used. 

ACKNOWLEDGEMENT 

This work was sponsored by the Stichting Nationale Computerfaciliteiten (National Com- 
puting Facilities Foundation, NCF) for the use of supercomputer facilities, with financial 
support from the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Netherlands 
Organization for Scientific Research, NWO). 

The authors wish to acknowledge the contribution of Mr. P. J. Hartc to a well-designed 
computer program for the calculation of the impurity resistivity. 



17 



REFERENCES 

1 J. Bass, in Landolt-Bornstein, Numerical Data and Functional Relationships in science 
and technology, New Series, edited by K.-H. Hellwege and J. L. Olsen (Springer- Verlag, 
Berlin, 1982), Vol. 15a. 

2 P. Dutta and P. M. Horn, Rev. Mod. Phys. 53, 497 (1981). 

3 P. M. Boerrigter, A. Lodder, and J. Molenaar, Phys. Stat. Sol. (b) 119, K91 (1983). 

4 N. Papanikolaou, N. Stefanou, and C. Papastaikoudis, Phys. Rev. B 49, 16117 (1994). 

5 1. Mertig, E. Mrosan, and P. Ziesche, in Multiple scattering theory of point defects in 
metals: electronic properties, edited by W. Ebeling, W. Meling, A. Uhlmann, and B. 
Wilhelmi (Teubner, Leipzig, 1987). 

6 T. Vojta, I. Mertig, and R. Zeller, Phys. Rev. B 46, 15761 (1992). 

7 J. van Ek and A. Lodder, J. Phys. Cond. Mat. 3, 7363 (1991). 

8 1. Mertig, R. Zeller, and P. H. Dederichs, Phys. Rev. B 47, 16178 (1993). 

9 J. van Ek and A. Lodder, J. Alloys Comp. 185, 207 (1992). 

J. P. Dekker, A. Lodder, and J. van Ek, Phys. Rev. B 56, 12167 (1997). 

1 J. M. Ziman, in Principles of the theory of solids, edited by J. M. Ziman (Cambridge 
University Press, Cambridge, 1972). 

2 A. Lodder and J. P. Dekker, in Proceedings of the First International Alloy Confer- 
ence(Athens, 1996), edited by A. Gonis, A. Meike, and P. E. A. Turchi (Plenum, New 
York, 1997), pp. 467-477. 

3 A. Lodder and P. J. Braspenning, Phys. Rev. B 49, 10215 (1994). 

4 For a review, see A. Lodder and J. P. Dekker, Phys. Rev. B 49, 10206 (1994). 

5 J. Molenaar, A. Lodder, and P. T. Coleridge, J. Phys. F: Met. Phys. 13, 839 (1983). 

18 



P. M. Oppeneer and A. Lodder, J. Phys. F: Met. Phys. 17, 1901 (1987). 

R. H. Lasseter and P. Soven, Phys. Rev. B 8, 2476 (1973). 

R. Schdpke and E. Mrosan, Phys. Stat. Sol. (b) 90, K95 (1978). 

J. Friedel, Nuovo Cimento Suppl. 7, 287 (1958). 

R. O. Simmons and R. W. Balluffi, Phys. Rev. 117, 62 (1960). 

Calculations of the electromigration wind force for this series of impurities in Ag were 
published previously in Ref. 10. 



19 



FIGURES 





impurity E 




o 
o 




FIG. 1. Definition of the electric field directions in the FCC structure. 
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FIG. 2. Definition of the electric field directions in the BCC structure. The directions 
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dicular to Em are equivalent. 
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FIG. 3. Impurity resistivity of 3d and 4sp atoms in Al. For the 3d metals constructed potentials 
are used with either one (filled circles) or two (filled squares) 4s electrons. Results obtained for 
the As 2 4p n atoms are also indicated by filled squares. Open circles correspond to resistivity values 
obtained with shifted potentials for Ca(4s 2 ), the 3d™ 4s 1 transition metal atoms and the 4p atoms. 
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FIG. 4. Vacancy pair with two different orientations with respect to the current. The geomet- 
rical cross-section is smaller, when the vacancy pair is aligned with the current. 



21 




FIG. 5. Calculated and measured resistivities of 5sp impurities in Ag. 
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FIG. 6. Calculated resistivities of pairs of a 5sp impurity and a vacancy in Ag. The resistivities 
the single impurities are given for comparison. 
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FIG. 7. Calculated resistivities in Ag of 5sp impurities, located next to a vacancy (initial 
position) and at the saddle point position. 




FIG. 8. Calculated and measured resistivities of 3d impurities in V. 
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FIG. 9. Calculated resistivities of 3d, Ad and 5d impurities, located next to a vacancy in V. 
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FIG. 10. Calculated resistivities of 3d, Ad and 5d impurities, located at the saddle point position 
in V. 
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